10 Observables and Propagators 



In the past lectures we have considered the properties of a number of field the- 
ories that describe physical systems isolated from the outside world. However, 
the only way to investigate the properties of a physical system is to interact with 
it. Thus, we must consider physical systems which somehow are coupled to their 
surroundings. We will have two situations in mind. In one case, we will look at 
the problem of the interaction of states of an isolated system, i.e., a scattering 
problem. In this case we prepare states (or wave packets) which are sufficiently 
far apart so that their mutual interactions can be neglected. The prototype 
is the scattering of particles off a target or each other in a particle accelerator 
experiment. Here the physical properties are encapsulated in a suitable set of 
cross sections. In the second case we will imagine that we want to understand 
that properties of a large system "from outside" , and we will consider the role of 
small external perturbations. Here we will develop a general approach, known 
as Linear Response Theory. We will couch our results in terms of a suitable set 
of susceptibilities. This is the typical situation of interest in an experiment in a 
condensed matter system. 

In both cases, all quantities of physical interest will be derived from a suitably 
defined Green function or propagator. Our task will be twofold. First we will 
determine the general expected properties of the propagators, in particular their 
analytic properties. Second we will see that their analytic properties largely 
determine the behavior of cross sections and susceptibilities. 

10.1 The Propagator in Electrodynamics 

In a Classical Field Theory, such as Classical Electrodynamics, we can inves- 
tigate the properties of the electromagnetic field by considering the effect of a 
set of well localized external sources. These can be electric charges or, more 
generally, some well defined distribution of electric currents. The result is fa- 
miliar to us: the external currents set up a radiation field which propagates in 
space, a propagating electromagnetic field. In Maxwell's electrodynamics these 
effects are described by Maxwell's equations, i.e., the equations of motion of 
the electromagnetic field in the presence of a current distribution j^(x) 



where we have assumed the Lorentz gauge condition (d^A 11 = 0), i.e., the wave 
equation. 

In Classical Electrodynamics the solutions to this equation is found by using 
of the Green Function G(x, x'), 



8 2 A ft (x) =f{x) 



(1) 



d 2 x G{x,x') = 5 i {x-x') 



(2) 



which satisfies the boundary condition that 



G(x, x') = if x < x„ 



(3) 



f 



This is the retarded Green Function, and it vanishes for events in the past (i.e., 
for x < x'q). 

The wave equations in the presence of a set of currents j^(x) is an inhomo- 
geneous partial differential equation (p. d. e. ). For times in the remote past 
(i.e., before any currents were present) there should be no electromagnetic field 
present. The choice of retarded boundary conditions guarantees that the system 
be causal. The solution to the inhomogeneous p.d.e. is, as usual, the sum of 
an arbitrary solution of the homogenous equation, A? n , which represents a pre- 
existing electromagnetic field, and a particular solution of the inhomogeneous 
field equation. We write the general solution in the form 

A"(x) = A$ n {x) + J dUG R {x x')f(x') (4) 

where A^ n (x) is a solution of the wave equation in free space (in the absence of 
sources), 

d 2 A? n (x) = Q (5) 

Thus, all we need to know is the Green Function. Notice that the choice of 
retarded boundary conditions insures that, for x$ < x' , A fJ/ (x) — A^ n (x), since 
Gr(x — x') = for x < x'q. 

Let us solve for the Green Function G(x,x'). This is most easily done by 
considering the Fourier transform of G(x, x') 

where p =J>^ and p ■ x = p n x n — p ■ x. It is easy to check that the Fourier 
transform G(p) should be given by 

G{P) = ~J 2 (7) 

But we have two problems now. First, G(p) has a singularity at p 2 = Pg—p 2 = 0, 
i.e., at the eigenfrequencies of the normal modes of the free electromagnetic field 
Po = i\p\- Thus the integral is ill-defined and some definition must be given 
about what to do with the singularity. Second, this G(x,x') does not satisfy 
(at least not in any obvious way) the boundary conditions. We will solve both 
problems simultaneously. Let us define the Retarded Green Function by 

G R (x, x') = 6(x - x' )G(x, x') (8) 

where 



Therefore, Gr(x, x 1 ) vanishes for x < x' and is equal to G(x — x') for x > x' . 
The step (Heaviside) function 9(a;) has the formal integral representation (a 
Fourier transform) 
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Q(x) = lim / duj (10) 
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where the integral is interpreted as a contour integral over the contours of Fig. 
1. Thus, when closing the contour as in the case x < 0, the contour does not 
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Figure 1: Contour in the complex plane that defines the function <d(x). 

contain the pole and the integral vanishes. Conversely, for x > 0, we close the 
contour on a large arc in the upper half plane and pick up a contribution from 
the enclosed pole equal to the residue, e~ ex which converges to 1 as e — > + . 
Notice that the integral on the large arc in the upper half plane converges to 
zero (for arcs with radius R — > oo) only if x > 0. 

We define the retarded Green function Gr(x-x') by the following expression, 

f d 4 p e~ l P ' ( x ~ x ') 

G R (x-x') = -i lim / jttuI r~~T2 ~ ( U ) 

e -o+J (2tt) 4 (po + ie) A -p 2 

which satisfies all of the requirements. We can also define the advanced Green 
function Ga(x — x') which vanishes in the future, but not in the past, 

G A (x-x') = x a -x' Q >0 (12) 

by changing the sign of e, i.e., 

f d 4 p e^iP ' ( x ~ x ') 

Ga(x - x) = -i lim / , — ^ (13) 

e -o+7 (2tt) 4 {po-ie) 2 -p 2 
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Figure 2: Contour in the complex plane that defines retarded Green function. 

Now the poles are on the upper half-plane and we get zero for x' — xq > but 
we pick up the poles for x' — xq < when we close on the lower half plane. 

There are still two other possible choices of contours, such as the one shown 
in figure 3. In this case, when we close on the contour 7 + in the upper-half 
plane, we pick up the pole at —\p \ and the resulting frequency integral is 



dpo e~ iv ' ( x ~ x ') ' {x - x ') + i\p\(x - x' ) 

2^ Pl-P 2 = 2(-\p\) 



while in the opposite case we get 

dp^e-iP' {x - x') & ip- (x - x') - i\p\(x 



Tj- 2ni pl-p 2 2\p\ 
Thus, this choice the contour yields the Green function 



x Q < x' (14) 



x > x' (15) 



G F (x -x') = -j {eK - -0) e-^l^o - -0) + ip- (x-x') 

+e(x - x' )e +i \P\( x '° - x °) - (*"*')} 

(16) 
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Figure 3: Contour in the complex plane that defines Feynman propagator, or 
time ordered Green function. 

which is known as the Feynman Green function and it bears a close formal 
resemblance to the mode expansions of free field theory. We will see below that 
there is indeed a reason for that. Notice that, in Maxwell's theory the zero 
component for the 4- momentum p is equal to \p\, and that the integration 
measure is Lorentz invariant. Also notice that this Feynman propagator or 
Green function propagates the positive frequency modes forward in time and the 
negative frequency modes backward in time. The alternative choice of contour 
simply yields the negative of Gf{x — x'). 

10.2 The Propagator in Non-relativistic Quantum Me- 
chanics 

In non-relativistic quantum mechanics, the motion of quantum states is governed 
by the Schrodinger Equation 

(ihdt - = (17) 

where H = H + V. Let V be some position and time-dependent potential 
which vanishes (very slowly) both in the remote past (t — > — oo) and in the 
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remote future (t — > +00). In this case, the eigenstates of the system are, in 
both limits, eigenstates of H - If V(x, t) varies slowly with time, the states of H 
evolve smoothly, or adiabatically. Thus, we are describing scattering processes 
between free particle states Let F{x't'\xt) denote the amplitude 

A time 



outgoing wave 




Figure 4: A scattering process. 



F{x't'\xt) = (x't'\xt) (18) 

We have already introduced this amplitude when we discussed the path integral 
picture of Quantum Mechanics. 

Let us suppose that at some time t the system has is in the state \tp(t)) = \xt). 
At some time t' > t, the state of the system is \ip, t') which is obtained by solving 
the Schrodingcr equation 

ihd t \il>) = H\i>) (19) 
where H is, in general, time-dependent. The formal solution of this equation is 

-- f dt"H(t") 
m'))=Te h Jt \^(t)) (20) 



G 



where T is the time ordering symbol, i.e., 
-i /"*' 

— / dt"H(t") 00 1 /_A" /■*' Z"*— 1 
T e h Jt = ~ { [y) J *i ' ' ' J ^(h) . . . H(t n ) 

(21) 

Thus, the amplitude F(x't'\xt) is 

F{x't'\xt) = (x't'\xt) = (x'\ijj{t')) (22) 
Hence, if the initial state \tp(t)} is \tp(t)} = \x) we get 



F(x't'\xt) = (x'\Te n Jt \ x ) (23) 

For the case of an arbitrary initial state \ip(t)) we get 



. dt"H(t") 

r.'\ ft. 



I 

(x'm')) = (x'\e h Jt ' \m) (24) 

Since the states {\x >} are complete, we have the completeness relation 



dx\x)(x\ (25) 



which allows us to write 



,+oc -i f dt"H{t") 

/ dx'(x'\e h Jt 

J — OO 



(x'|#'))=| dx'(x'|e \x)(x\4>(t)) (26) 

so, we get 



dx'{x't'\xt)ip(x,t) (27) 

-oo 

In other words, the amplitude F(x't',xt) is the kernel of the time evolution for 
arbitrary states. The amplitude F{x't'\xt) is known as the Schwinger function. 

The initial state \ip(t)} and the final state \ip(t')} are connected by the evo- 
lution operator U (f, t) 

- l - f dt"H(t") 
U(t',t) = Te n Jt (28) 

which is unitary since, as a result of the Hermiticity of the Hamiltonian, we have 

i 



J dt"H(t") 



U\t',t) = Te n Jt 

- l - fdt"H{t") 
= Te h Jt> =U(t,t') 



(29) 
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Be definition U(t,t') is the inverse of U(t',t) since it involves the states back- 
wards in time. In addition, the operator U(t',t) obeys the initial condition 

lim £/(i',*) = J (30) 

where I is the identity operator. If the Hamiltonian is time-independent, the 
evolution operator is 



-1 f dt"H{t") >:/' 
U{t',t)=Te h Jt = e h V ' (31) 



These ideas will also allow us to introduce the Scattering Matrix (or S -Matrix). 
If ipi(x,t) is some initial state and ipf(x',t') is some final state, the matrix 
elements of the ^-matrix between states tpi and ipf,Sfi are obtained by evolving 
the state fa up to time t' and projecting it onto the state ipf. Namely 

S fi = lim lim [ dx [ dx'i>f(x't'){x't'\xt)i>i(x,t) (32) 

Let us define the Green function or propagator G(x't'\xt) 

G(x't'\xt) = - l -Q{t' - t)(x't'\xt) = ^Q{t' - t)F(xt'\xt) (33) 

It satisfies the equation 

(ihd t > - H(x')) G(x't'\xt) = 6(x - x')S(t - t') (34) 

with the boundary condition 

G{x't'\xt) = if t' < t (35) 

Hence G(x't'\xt) is the retarded propagator. 

In terms of G, the S'-matrix is given by (recall that t' > t) 

S fi = i lim lim / dx / dx'ip* f (x' , t')G(x't'\xt)ipi(x, t) (36) 

t'— t+oo i^- oo J J 

Let us consider now the case of a free particle with Hamiltonian Hq which is 
coupled to an external perturbation represented by a potential V(x, t). The free 
Green function, Go, satisfies the equation 

(ihd t > - H (x')) G = S(x' - x)5{t' - t) (37) 

Go can be regarded as the matrix elements of the following operator 

G (x't'\xt) = (x't'\(ihd t - Hy^xt > (38) 

Clearly, G satisfies the same equation but with the full H, i.e., 

(ihd t > -H{x'))G = 1 (39) 
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Hence, we can write 

[(ihd t > - H ) - V] G = 1 
By using the definition of Go, we get the operator equation 

(Go" 1 - V)G = 1 

Thus, G satisfies the integral equation 



(40) 
(41) 



G(x't',xt) = G {x't'\xt) + / dx" / dt"Go(x't'\x't")V(x't")G(x't"\xt) (42) 



which is known as the Dyson Equation. It has the formal operator solution 



G _1 = Gn 1 - V 



(43) 



The integral equation can also be solved by iteration, which amounts to a per- 
turbative expansion expansion in powers of V. The result is the Born series. 
Using an obvious matrix notation we get 



G = Go + G a VG a + G a VG VG + . . 



(44) 



We will represent this series by a set of diagrams. Let us consider the first term 
to which we assign the diagram of figure 5. The oriented arrow ranging from 



A 



time 



(x',t>) 



(x,t) 



space 



Figure 5: The zeroth order term. 
(xt) to (x't') represents the unperturbed propagator Go(x't'\xt). Because of 
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the causal boundary conditions obeyed by Go it can only propagate forward in 
time. The second term, of the series, the Born approximation, SG^ 

5G {1) (x't'\xt) = J dx' J dt"G (x't'\x"t")V{x"t")G (x"t"\xt) (45) 

is represented by the diagram where the "blob" represents the action of the 

t 

(x',t>) 




{x",t" 



Figure 6: The first order term: the Born Approximation. 

potential V. In general we get a diagram of the form which represents a multiple 
scattering process. Notice that all the contributions propagate strictly forward 
in time. 

Let us compute the propagator Go(x't'\xt) for a free spinless particle in 
three-dimensional space. The Hamiltonian H is just. H = — ^V~- Thus. Go 
obeys the equation 



{%m> + ^-vl)Go(x't'\xt) = 6 3 (x' - x)S(t' - t) 



(46) 



we causal boundary conditions, G(x't'\xt) = (if t' < t). Given the symmetries 
of this very simple system we can Fourier expand G . 

G (x't'\xt) = J^J ^Go&^ei^'-^-^'-V (47) 



By direct substitution we find that Go(p, ui) must be given by 

1 



G (p>) 



(48) 



2m 
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space 



Figure 7: Multiple scattering processes. 



p 2 

Notice that, once again, Go(p, w) has a pole at ftw — which is the dispersion 

2m 

law (i.e., , the "mass shell"). Apart from being singular, this "solution" does not 
obey the causal boundary condition. We will enforce the boundary condition 
by deforming the integration contour in the complex frequency plane. Following 
our previous discussion on the Green function for classical electrodynamics we 
simply move the pole by an infinitesimal positive amount e into the upper half- 
complex frequency plane. We write the retarded propagator as 

GT(P» = 1~ 2 (49) 

hw - h ie 

2m 

and we will take the limit e — > + at the end of our calculations. Thus the 
frequency integral is equal to 



tiuj — h ie 

2m 

Hence the full Green function Go(x't'\xt) is 



i 



d 3 P A£=v 



G {x't'\xt) = --@{t' -t) J h2m y h r (51) 
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Let ipp(xt) denote the wave functions for the stationary states \p), i.e., 

1 -tP ■ x+ T E (p)t 
M^)^^-^e h P ft (52) 

P 2 

where E(p) = - — . We see that Go(x't'\xt) can be written in the form 
2m 

G (x't'\xt) =- % - Q(t' -t) j d 3 p ^{x't')^{xt) (53) 

In general, if the Hamiltonian has a complete set of stationary states {\n)} with 
wave functions ip n (x) and eigenvalues E n , the Green function is 

(y°\x't'\xt) = - l -e(t'-t)J2Mx't'W n (xt) (54) 

n 

where 

— "t E nt 

ipn(xt) = ip n {x)e n ' (55) 

If the system is isolated, the Hamiltonian is time-independent and the Green 
function is a function of t — t'. In this case it is convenient to consider the 
Fourier transform 

Gl et (x';x\t' -t)= I" —G™ t (x',x;u J )e lLU ( t - t 'y h (56) 

J-oo 27r 

where we have to pick the correct integration contour so that Gjf* is retarded. 
Quite explicitly we find 

(*'; x\u>) = \ lim V tll^fiM (57) 

■a n 

The denominators in this equation have zeros on the real frequency axis as 
e — > + . Thus, the Green function is a series of distributions. In the limit 
e — > + we can write 

lim — !— =V-- in5(x) (58) 
e^o+ x + ie x 

where where V denotes the principal value and S(x) is the Dirac delta function. 
Hence we can write the following expressions for the real and imaginary parts 
of the Green function 

n n 

lmGl ct (x';x\uj) = 

n 

(59) 
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We can use these results to find an expression for the density of states p(ui) 

p (uj) = ]T 5(hw - E n ) (60) 

n 

in terms of the retarded Green function of the form 

p(uj) = ~lm J dx G T et (x;x\u) (61) 

In other words, the spectral density is determined by the imaginary part of the 
Green function. Below we will find similar relationships for other quantities of 
physical interest. 

Let us close by noting that there is a close connection between the Green 
function in Quantum Mechanics and the kernel of the diffusion equation. In 
d-dimcnsional space the Green function is 



G (x',t'\x,t) = - l -Q(t'-t) J 



(27r) d v ; 



By completing squares inside the exponent 



2m 



{t' -t)-$ -{x" -x) 



(t'-t) 
2m 

t'-t 



_ x — x 



2m 



x — x 
t'-t 



in 



mix' — x\ 2 



2\t'-t\ 



(63) 



we can write 

G (x',t'\x, t) 



d d p 



t'-t 
2m 



t'-t 



in 



2h 



rn- 



t'-t 



(64) 



After a straightforward integration, we find that G is equal to 



G (x',t'\x,t) = --e(t' -t) 



mh Y /2 c 2h (t'-t) 

2n(t' -t)J 



im \x — x\ 



(65) 



This formula is strongly reminiscent of the Kernel for the Heat Equation (or 
Diffusion Equation) 

drip = D\j 2 iP (66) 

where D is the diffusion constant and r is the diffusion time. Indeed, after 
an analytic continuation to imaginary time, t — ► ir, the Schrodinger equation 
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becomes a diffusion equation with a diffusion constant D = The Green 

function (or Heat Kernel) for the Heat equation is 

\x' - x\ 2 

K(x'r'\xr) = 6(r' - r)-^-^ e ~ ">(t> - t) (6?) 

which, of course, agrees with Go- This connection between quantum mechan- 
ics and diffusion processes is central to the path-integral picture of quantum 
mechanics. 

10.2.1 Analytic Properties of the Propagator of the Relativistic Real 
Scalar Field and of the Dirac Field 

The propagator for a free real relativistic field Gq(x,x') — — i(0\T(j)(x)(j)(x / )\0} 
is the solution of the partial differential equation 

(d 2 +m 2 ) G Q (x,x') = 5 4 (x-x') (68) 

which we have discussed and solved before. Go(x,x') can be calculated by the 
usual Fourier expansion methods, i.e., 

Go(x,x') = / Go(p) e +i *>< x - x ') (69) 

where Go{p) is given by 

G (p) = -^-5 (70) 

where p 2 = p^p^- Once again, we will have to give a prescription for going 
around the poles of Go{p) which yields the correct boundary conditions. The 
poles of Go{p) are located at po = ±\/p 2 + m 2 ^fie. We will use the integration 
paths shown in figure 8. On the integration path 7+ we find the result 



dpo e +i Po( x o - x 'o) _ e -iy/f 2 ~Tm 2 (x - x' ) 

+ 2tt p 2 - (p 2 + m 2 ) + ie e ^o+' % 2^p 2 + m 2 ^ 



provided xq > x' Q . Similarly, on 7 we get 



a ip (x - x' ) JVp 2 +m 2 (x - xq) 



dp e l P^ x °- x o) e % 



' 7 _ 2tt p 2 - (p 2 +m 2 ) +ie e ^o+* % 2y/p 2 ' ™ 2 



provided xo < x' Q . By collecting terms we get 
G (x,x') = *6(x -x ) / { 2^L{p) e+ ^' { *~*' ] X '° ] 



(73) 



14 



Im po 




Figure 8: Wick Rotation. 



This formula shows that Gq{x, x') does satisfy the required boundary condition. 
It also shows that the positive frequency components of the field propagate 
forwards in time while the negative frequency components propagate backwards 
in time. 

The simplest way to compute Go(x, x') is by means of an analytic continua- 
tion (or Wick rotation) to imaginary time, xo — > ix^. This amounts to a rotation 
of the integration contour from the real po axis to the imaginary po axis, namely 
Po — » ipi- The Wick-rotated or Euclidean Green function Gq(x,x'), i.e., the 
Euclidean propagator that we calculated before, is given by 



d A p e _<p ( Xo_x o) 



= * / ,oJ>4 2^ 2 = iG*{x,x') (74) 



where p 2 = - Yh=i V%V% and p ■ x = - Y.t=i Vi%i- 

The time-ordered, or Feynman, propagator does not obey causality since it 
does not vanish for space- like separated events (s 2 < 0) . We can define a causal, 
or retarded, propagator which obeys the causal boundary condition G(x, x') = 
except inside the forward light-cone. Similarly, it is also possible to define an 
advanced, or anti-causal, propagator which vanishes outside the backward light- 
cone. We will discuss only the retarded propagator G let (x — x'). 

The retarded propagator is defined by computing the frequency integral on 
the path shown in the figure 9, G rc t(a; — x') is given by 
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Figure 9: Integration path for the retarded Green function. 



r rfin ■- <-> ->\ r + °° rinr, piPo( x O - x' ) 

G iet (x-x') = - AP e -*P-(x-x) / |^ ° 

J (2tt) 3 J_ 00 2tt p^-u> 2 (p) 

/■ rf 3 p -jg- (f - x 1 ) ^ojp (x - 4) 



Tf2oj(p) 



p -uj{p)-ie p + Lu(p)-ie 



(75) 



Hence, we get 

G Ict (x - x') = -iO(x Q - x' ) 

f d3p \ Joj(p)(x - x' ) - ip- (x-x) _ -iu>(p)(x - x' ) + ip- (x - f ) 
J (2n)*2u(p) I 

(76) 



The integral over the momentum variables is just the quantity iA(x — x'), that 
we have encountered before, 

t A(x - x') = Wx), <f>(x')] ee (0| [<j>(x), cj>{x')] |0) (77) 

which allows us to write 

G rct (x - x') = e(x - x' )A(x - x') = -iQ{xo - ^)(0| [<t>(x), <t>{x')} |0) (78) 

The same line of argument we have used here for the scalar field can be used 
for the Dirac field. The vacuum state of the Dirac theory is defined by filling 
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up all negative energy states. We now imagine that an external potential is 
adiabatically switched on and that, before that happened, there was an electron 
e— propagating in the system. If the potential is not too strong, we can still 
describe its effects by means of a Born series of multiple scattering processes 
(Fig.10). 




Figure 10: Scattering processes of a Dirac field. 

The electron scatters off the potential at 1 and as a result it propagates up 
to 3. If the potential has the correct matrix elements, at 3 the positive energy 
solution may turn into a negative energy state. In general this won't be allowed 
since all negative energy states are filled up, unless a negative energy state got 
emptied (by the action of the potential) before the electron became scattered 
into that state. Indeed, the potential can create a pair out of the vacuum, e.g., 
Fig.ll: 



A 




space 



Figure 11: Pair creation. 

Thus, it can remove an electron from an occupied negative energy state and 
to promote it into a previously empty positive energy state. Hence, if the time t 
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lies between i 2 and £3, there are three different states propagating in the system 
(Fig.12): 

A 




space 



Figure 12: An intermediate state with an electron-positron pair. 

1. a positive energy state which disappears at 3 

2. a negative energy state which propagates backwards in time from 3 to 2 

3. a positive energy state which appeared at 2. 

An alternative interpretation is that an electron-positron pair was created at 2 
and that the positron annihilated with the original electron at 3. This process 
clearly shows that the Dirac theory is a quantum field theory, and cannot be 
described within the framework of quantum mechanics with a fixed number of 
particles, as in the non-relativistic case. Thus, the Fock space description is 
essential to the relativistic case. 

These arguments suggest that we may want to seek a propagator which 
propagates positive energy states forward in time while negative energy states 
propagate backwards in time. This is the Feynman, or time-ordered, propagator 
Sp a (x,x r ). It is straightforward to see that these requirements are met by 

S F a \x,x') = -i(0\Tip a (x)ij a ,(x')\0) (79) 

which satisfies the equation of motion 

{i$-m)S F {x,x') = 5 i {x-x') (80) 

The same methods that we used for the scalar field yield the answer (dropping 
the spinor indices) 

S F (x, x') =S F (x- x') = f J + ™ . e-*P ■ ( x - x ') (81) 

v ; y 'J (2tt) 4 p 2 - m 2 + it K ' 
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where an ie has been introduced in order to get the correct boundary conditions. 
We have also shown that S F (x — x') satisfies 

S F {x-x') = -{i$ + m) / j^— Y — = {i$ + m)G F {x-x l ) (82) 

J (2-7T) 4 p — 772 -|- ie 

where G F (x — x') is the Feynman or time-ordered propagator for a scalar field. 
The ie prescription insures that positive energy states propagate forward in time 
while negative energy states propagate backward in time. 



10.2.2 The Propagator for the Non-Relativistic Electron Gas 

Let us now discuss the propagator, or one-particle Green function, for a non- 
relativistic free electron gas at finite density and zero temperature. It is defined 
in the usual way 

G% a '(x,x') = -i<gnd|T^ a (z)^,(z')|gnd) (83) 

where aa' are spin indices. This propagator can be used to compute a number of 
quantities of physical interest. For example, the average electron density (n(x)) 
is 

(n(x)) = (i>*(x)i>(x)) = -i Jim_ lim tr G(x,x') (84) 
Likewise the current (j(x)} is 

n 



(M) = ^rtr (^(x) (v^(i)) - (y x ^{x))^{x)) 

(vj-V f -)trG(x,a;') 



= lim lim 

771 x' — »x x' ^x a 



G(x,x>) = J^J^Gip^e^-V-^-^-t')] 
G(p) e ip - (x - x "> 



I 



d 4 p 



(85) 



and the magnetization M(x) = {tp' { a {x)a a ^ p(x)) , 

M(x) = -i ]m lim tr [G(x, x)a] (86) 

x' — >x x' D — >XQ 

Let us compute G(x,x') by the standard method of Fourier transforms. 

d 3 p 

(27T) 



(87) 



where po = u> and t = x . The field ip(x) can also be expanded (h= 1) 

m = / m ^ ^-*- Em (88) 
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since ip(x) satisfies the equations of motion. Here E(p) are the single particle 

2 

energies, E(p) = ^ — A*. Thus, we get 

- -z/^/^eW(gnd W p)^(q)|gnd) e i (P---^1*) 

+ (03©(-i)<gnd|V t (q)^(p)|gnd)e i ^-^-^) 

(89) 

Recall that |gnd) is the state in which all negative energy states (E(p) < or 
e(p) < n) arc filled. Thus 

G(x,t) = -i J ^[QWil-nW-Qi-tfrme'V-*- 15 ^ (90) 
where n(p) is the (zero-temperature) Fermi function 

"« = { J <w 

Hence, we can write the Fourier transform G(p, u) in the form 

G(p,w) = -ije(|£|-pjr)jT e i{ - u - E< ^*dt-Q{pF-\p[) j cfe-^- B(p))t J 

(92) 

The integrals in this expression define distributions of the form 
dte^ = lim / dte ist ~ et 



= i lim — = i \V ine(s) 

e^0+ s + iA \ s 



(93) 



where V- is the principal value 



V- = lim -^i-. (94) 
s e^a s 2 + e 2 

and <5(s) is the Dirac 5- function. 

We can use these results to write an expression for G (p,uj) 

dUM.-aga-eU-^dSL (95) 

U! - E(p) + 16 UJ - E(p) - IE 

where E{p) = e(p) — \i. An equivalent (and more compact) expression is 

G Q (p,w) = , 1 . — — (96) 

u - E(p) + le sign(|p| - p F ) 
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Notice that 



ImG (p» = -7t6(|p1 -p F )5(u> - e(p) + p) +7r9(p F - \p\)5(w - e(p) + p) 
= -nS(u; - e(p) + p) [9 (\p\ - p F ) - Q( PF - \p\)] 

(97) 

The last identity shows that 

sign ImGo(p, w) = —sign w (98) 
Hence, we may also write Go (p, w) as 

1 



G (p,uj) 



lu — E(p) + ie sign lu 



(99) 



With this expression, we see that Go{p,w) has poles at u = E(p) and that all 
the poles with iv > are infinitcsimally shifted downwards to the lower half- 
plane, while the others are raised upwards to the upper half-plane by the same 
amount. Since E(p) = e(p) — /x, all poles with e(p) > n are shifted downwards, 
while all poles with e(p) < p, are shifted upwards. 



Im i 



Re U 



\\x\\\\\x\\x\\\i4\x\\x\\x.\.\\.\.\\.\.\\\.\\\.\.\ 
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Figure 13: Analytic structure of the Green function at finite chemical potential 
p. 
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10.3 The Scattering Matrix 

The problems of real physical interest very rarely involve free fields. In general 
we have to deal with interacting fields. For the sake of definiteness we will con- 
sider a scalar field but the ideas that we will discuss have general applicability. 

10.3.1 Interactions and the Interaction Representation 

The field <p(x, t) in the Heisenberg representation is related to the Schrodinger 
operator <j>{x) through the time evolution generated by the equations of motion, 
i.e., 

<f>(x, t) = <f>(x) = e iHt (j){x, 0) e- lHt (100) 

(for h = 1). Recall that in the Schrodinger representation the fields are fixed 
but the states evolve according to the Schrodinger equation 

H\$) = id t \<f>) (101) 

Conversely, in the Heisenberg representation, the states are fixed but the fields 
evolve by following the equations of motion, i.e., 

id t <t>(x) = [<j)(x),H] (102) 

In an interacting system, the problem is precisely how to determine the evolution 
operator e lHt . Thus, the Heisenberg representation cannot be constructed a 
priori. 

Let us assume that the Hamiltonian can be split into a sum of two terms 

H = H + H int (103) 

where Hq represents a system whose states are fully known to us (i.e., a problem 
that we know how to solve) and H int represents the interactions. For technical 
reasons we will have to assume that Hi nt (t), as a function of time, vanishes (very 
smoothly) both in the remote past and in the remote future. We now define 
the Interaction Representation. In the Interaction Representation, the fields 
evolve as free Heisenberg fields, i.e., 

id t cj ) - m {x) = [<i) in ,H ] (104) 

These operators create and destroy free incoming states (incoming since as t — > 
— oo there are no interactions). In the absence of interactions, the states would 
not evolve but, if interactions are present, they will. There is a unitary operator 
U (t) which governs the time evolution of the states and the S'-matrix. 
We want to find an operator U (t) such that 

4>(x,t) = U- 1 <j) in {x,t)U{t) (105) 

where <j>(x,t) is the Heisenberg operator. The operator U(t) must be unitary 
and satisfy the initial condition 

lim U{t) = I (106) 

t — > — OO 
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Since U(t) is unitary and invertible, it must satisfy the condition U 1 (t) = 
W{t). Thus 

±U{t)U-\t) + U{t)jU-\t) = Q (107) 
Since the operator 0; n obeys its equation of motion, we get 

d t <f) in = d t U tftU^ + U d t <j> U^ + Ucp dtU^ 1 

= d t U(l)U- 1 + iU[H,<j>}U- 1 + U(j>U- 1 UdtU- 1 

(108) 

In other terms, 

dtfrn = iU(t)[H(<f>), <P]U-\t) + dtUU- 1 ^ + frnUdtU- 1 (109) 

Similarly, the Hamiltonian must obey the identity 

H(cj> iTl ) = U(t)H((t>)U- 1 (t) (110) 

which implies that should obey 

dtcpin = i[H(<f> in ),(j> in } + [(d t U) U-\ fa] (111) 

Since </>;„ obeys the free equation of motion, we find that the evolution operator 
U (t) must satisfy the condition 

i[#int(<&n), 0in] + [8 t U U~\ fa] = (112) 

Hence 

[i#int(0i») + (d t U) U~\ fa] = (113) 

for all operators <j)- m . Therefore this operator must be a c-number {i.e., propor- 
tional to the identity). Since lim^-oo -f/i n t(0in) = and \im t ^-ooU(t) = I, 
this c-number must be equal to zero. We thus arrive to an operator equation 
for U(t) 

id t U = H int (<k a ) U(t) (114) 

The operator U governs the time evolution of the states in the Interaction 
Representation, since the state |$)i n becomes 

[/(t)|$) in = |$(t)) (115) 

In particular, the outgoing states |$) ut (i.e., the states at t — > +oo) are also 
free states which are related to the in-states by the operator U (t) in the limit 
t — > +oo 

|$) out = lim tf(t)|$>i» = S\$) in (116) 
where S = lim t ^ +00 U (t) is the S -matrix. 
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The equation of motion for U(t) can formally be solved using methods that 
we have discussed before. The solution is 



— i 



dt'H int (<f> in (t')) 



U(t)=Te (117) 

where T is the time-ordering operator. In terms of the interaction part £- mt of 
the Lagrangian we get 

/+oo „ 
dt'H int (<p in (t')) i / d 4 x£ int (fr n ) 
,,,,, >„,-., -oo =T e J (118) 

10.3.2 Physical Information Contained in the S'-Matrix 

Let us compute transition matrix elements between arbitrary in and out states. 
Let \i, in) be the initial incoming state and \ f, out) the final outgoing state. The 
transition probability Wi^f is then given by 

W^ f = |(/,out|z,in)| 2 = |(/,in|,S|z,in)| 2 (119) 

since 

(/,out|5=(/,in| (120) 
We can split S into non-interacting and interacting parts 

§ = I + if (121) 

where / (the non-interacting part) is the identity operator, and T, the T-matrix 
represents the interactions (not to be confused with the time-ordering symbol!). 
In terms of the T-matrix, the transition probability is 

w^ f = \(m +l (f\mf (122) 

Although from now on we will discuss the case of a scalar field, the arguments 
can be generalized to all other problems of interest with only minor modifica- 
tions. 

Let us consider the situation in which the initial state |z, in) consists of two 
wave packets which only contain positive frequency components 

|2 ' in> = /^F / 2^(£)3^^)^^)^^ 2 ' in ) ( 123 ) 

The incoming flux is equal to J 2^(2^)3 \ f(p)\ 2 - Each component \pi,p2;in) will 
have a matrix element with the final state |/;out). Since we have translation 
invariance, the total 4-momentum should be conserved. If we denote by P/ the 
momentum of the state |/, out), we can write 

(.f\T\pip 2 ) = (2n) i 5 4 (Pf-pi-P2)(f\T\ Pl p 2 ) (124) 
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where T is called the reduced operator. T acts only on the energy shell, p 2 = 
m 2 . If we neglect the forward scattering contribution, the transition probability 
Wi^f becomes 

f d 3 Pl f d?p 2 f Ai [ d 3 q 2 
~* 7 2p?(2tt)3 7 2^(2tt)3 7 2g»(27r)3 7 2^(2tt)3 

(2n) 4 S 4 ( Pl +p2- qi - eft) (2tt) 4 <5 4 (P / - Pl - pa) 
/r(Pi)/2(P2)/i(gi)/ 2 (92)(/|T|p 1 p 2 )*(/|T| 9ig2 ) 

(125) 

The incoming states are assumed to be sharply peaked around some momenta 
Pi,p 2 , with a spread Ap, such that we can approximate 

(f\T\ PlP2 ) « (/|T| gi (ft) « (/|T|pip2) (126) 

Under these assumptions, the form for the transition probability is 

WW = / rf 4 x|/ 1 ( a; )| 2 |/ 2 ( a ;)| 2 (2^) 4 ( 5 4 (P / -pi -p 2 )|(/|T|p 1 p 2 )| 2 (127) 

where /(#) is the Fourier transform of /(p). The integrand of Eq.(127) is the 
transition probability per unit time and volume 

^g- = \f 1 (x)\ 2 \f(x)\ 2 (2n) 4 5\P f - Pl -p 2 )\(f\T\ Pl p 2 )\ 2 (128) 

In position space, the flux is i J d 3 xf* (x)dof(x). If f(x) is sufficiently smooth, 
we can approximate 

if*(x)d Q f(x)*2p \f(x)\ 2 (129) 

Let us assume that particle 1 is incident in the laboratory and that particle 2 
is at rest in the lab. The density of particles in the target is 

% = 2p8l/2(*)| a (130) 

where p 2 = m 2 since particle 2 is at rest. The incident flux is the velocity times 
the density of particles in the beam. Hence the incident flux $ in is 

$i„ = |^2p?|/ 1 ( a; )| 2 = 2^11/! (x)\ 2 (131) 

The differential cross section da is related to the transition probability by the 
relation 

dW^f dn 2 



dtdV dV 

Hence 



• $ in • da (132) 



2m 2 \f 2 (x)\ 2 2\p 1 \\f 1 (x)\ 2 da = 

\f 1 (x)\ 2 \f 2 (x)\ 2 (2n) 4 S i (P f -p 1 -p 2 )\(f\T\p 1 p 2 )\ 2 

(133) 
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Therefore, the differential cross section da is 

da = {2k)H\P } -p 1 -p 2 )-±-^f\T\p 1 p 2 )\ 2 (134) 

The quantity in the denominator TO2I.P1I can be written in the relativistic in- 
variant way 



m 2 \p\ =m 2 \jpf -m\ = [(p 2 • pi) 2 - m\m 2 2 ] 1/2 (135) 

Thus far we have not made any assumptions about the nature of the final state. 
If the process that we consider involve two particles going into n particles, the 

2 Scorning n outgoing 

particles particles 



Figure 14: A2^n scattering process. 

total differential cross-section becomes 

1 f d 3 p 3 f d 3 p n 



da _ 1 f d 3 p 3 f d A p n+ 2 

4 [(pa - mjmlf 2 J a (27r) 3 2p° ' ' 7a (2n) 3 2p° n+2 

n+2 

x \{p 3 , . . . ,p n+2 \T\p lP2 }\ 2 (2ir) 4 6(pi +p 2 - ^Po) 



(136) 



where A is an energy-momentum resolution. This expression shows that the 
central issue is to compute matrix elements of T. 
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10.4 Asymptotic States and the Analytic Properties of the Propa- 
gator 

We will show now that the S'-matrix elements can be calculated if we know the 
v.e.v. of time-ordered products of field operators (j){x) in the Heisenberg repre- 
sentation. This is the LSZ (Lchmann, Symanzik and Zimmermann) approach. 

In order to make this connection, it is necessary to relate the interacting 
fields <p(x) with a set of fields which create (or destroy) the actual states in the 
spectrum of the system, when these particles are far from each other. In any 
scattering experiment, the initial states are sharply peaked wave packets which 
can be constructed to be arbitrarily close to the eigenstates of the system. Of 
course, the true eigenstates are plane waves and any two such states will have a 
non-zero overlap in space. But wave packets which are essentially made of just 
one state will not overlap if the wave packets (the "particles" ) are far apart from 
each other at the initial time. Since they do not overlap, they do not interact. 
In this sense the spectrum of incoming states can be generated by a set of free 
fields. We will make the further assumption that the in and out states, the so- 
called asymptotic states, are created by a set of free fields. We will denote them 
by </>in(aO, and have the same form as the fields that appear in the Lagrangian. 
If the operators create the actual eigenstates, (f>i n {x) must be a free field. 

These assumptions amount to say that the states of the fully interacting 
theory are in one-to-one correspondence with the states of the non-interacting 
theory. In some loose sense, this hypothesis implies that the information that 
we can obtain from perturbation theory is qualitatively correct. 

However, these assumptions can fail in several possible ways. A mild failure 
would be the appearance of bound states which, of course, are not present in 
the unperturbed theory. This situation is actually rather common and it can be 
remedied without too much difficulty. An example of this case is the Landau 
theory of the Fermi liquid and its collective modes. There are however at least 
two ways in which this picture can fail in a rather serious way. One case is 
the situation in which the fields of the Lagrangian do not describe any of the 
asymptotic states of the theory. An example of this case is Quantum Chro- 
modynamics (QCD) whose Lagrangian describes the dynamics of quarks and 
gluons which are not present in the asymptotic states: quarks are confined and 
gluons are screened. The asymptotic states of QCD are mesons, baryons, and 
glue-balls which are bound states of quarks and gluons. 

Another possible failure of this hypothesis is the case in which the states 
created by the fields of the Lagrangian are not the true elementary excitations 
but rather they behave like some effective composite object of some more el- 
ementary states. In such case the true one-particle states may be orthogonal 
to the states created by the fields of the Lagrangian. This is a rather common 
situation in theories in 1 + 1 dimensions whose spectrum is generated by a set 
of soliton-like states, which are extended objects in terms of the bare fields of 
the Lagrangian. Something very similar happens in the theory of the fractional 
quantum Hall states. 

In what follows, we will not consider these very interesting situations and 
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assume instead that these hypotheses (or scenarios) actually hold. For the sake 
of concreteness, we will deal with scalar fields, and we will require to satisfy 
the free Klein-Gordon equation 



(d 2 + m 2 ) & n = 



(137) 



where m is the physical mass. In general, the physical mass is different from the 
mass parameter mo which enters in the Lagrangian of the interacting theory. 
As time goes by, this initial state evolves and the particles approach each other 
and begin to overlap. Interactions take place and, after some time, the system 
evolves to some final state which consists of a set of well defined particles, the 
out state. The unitary operator that connects in and out states is precisely 
the S'-matrix of the interaction representation. The only difference here, resides 
in the fact that the in and out states are not eigenstates of some unperturbed 
system, but the actual eigenstates of the full theory. 

This picture assumes that there is a stable vacuum state |0) such that the 
observed particles are the elementary excitations of this vacuum. The free fields 
4>i n (x) are just a device to generate the spectrum and have no real connection 
with actual dynamics. On the other hand, the interacting field <p(x) creates 
not only one-particle states but also many particle states. This is so because its 
equations of motion are non- linear. Hence the matrix element of 4>(x) and (j)i n (x) 
between the vacuum |0) and in-one-particle states |1), are generally different 
since </>; n creates only one-particle states. This must be true even as t — > — oo. 
We state this difference by writing 



The proportionality constant Z is known as the wave-function renormalization. 
If Z =/= 1, the operator 4>{x) must have a non-zero multi-particle projection. 
Notice that this is an identity of these matrix elements only. It is not an identity 
between the fields themselves. 

In the interaction representation it is possible to derive a similar looking iden- 
tity which originates from the fact that the unperturbed and perturbed states 
do not have the same normalization. It is important to stress that this approach 
makes the essential assumption that the states that are reached through pertur- 
bation theory in the interaction representation can approximate with arbitrary 
precision all of the exact states of the theory. This assumption is the hypothesis 
that the asymptotic states are generated by free fields. 

The physical asymptotic states satisfy canonical commutation relations and 
the commutator of a pair of much fields is 



where m is the physical mass. On the other hand, the interacting fields satisfy 




(138) 



(Opin(a;), <M*')]|0> = *A(x - x';m) 



(139) 



(0|[^),^')]|0) = ]T [(O|0(x)|n)(n|0(x')|O) - (; 



x <-» x')] 



(140) 



n 
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where {\n)} is a complete set of physical (in) states. The operators 4>{x) are 
related to the operator <fi(0) (i.e., at the origin at some time Xo = 0) by 

4>(x) = e ip - x c/)(0)e- ip - x (141) 

where P^ is the total 4- momentum operator. If we let P% be the 4- momentum 
of the state \n), we can write 

(0|[#r),</>(z')]|0) = £ [(0\m\n)e- w -< x - x ')<n|0(O)|O> - (a; - .*')] (142) 
We now insert the identity 

1 = J d 4 QS 4 (Q-P n ) (143) 

to get 

(0\[<P(x),cj>(x')]\0) = / d 4 Q ^ 4 (Q-P„)|<0|«K0)|n>| 2 ( e -«H*-*') - e «M*-*')) 

(144) 

We can summarize this result in terms of a spectral density p(Q) 

(0\[4>(x)^(x')]\0) = J -j^P(Q)(e- iQ<X - X 'W^"*')) (145) 
where p(Q) given by 

p(Q) = (2^) 3 SHQ P n )\(0\m\n)\ 2 (146) 

n 

Let us recall that A(x — x'\m) is given by 

= j^(Q a )S 4 (Q 2 - m 2 )e-^-^ 

(147) 

where e(Q°) = sign(Q°). Thus we can write 

(0Mx)^(x')]\0) = J ^ p (Q)e(Q°)e-^-^ (148) 

Since p(Q) is Lorentz invariant by construction, it can only be a positive function 
of Q 2 

p(Q) = <j(Q 2 ) > (149) 

Hence 

(O|^(x),0(x')]|O) = J ^a(Q 2 )e(Q°)e-^-^ (150) 
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We will now rewrite this expression in the form of an integral over the spectrum. 
Let us insert the identity 



to obtain 

(0|[^),^')]|0) 



1= / dfi 2 S(Q 2 - fi 2 ) 
Jo 



(151) 



<y(Q 2 HQo)e 



-iQ-{x — x') 



where a(fi 2 ) is 



(152) 
(153) 



a^ 2 ) = (2nrJ2S(P 2 n -^)\(0\m\n)\ 2 

n 

Thus, we can write the v. e. v. of the commutator as 

(0\[<p(x),(f)(x')]\0} =i / dn 2 <j(n 2 )A(x-x';m' 2 ) (154) 
Jo 

If we assume that the theory has a physical particle with mass to (i.e., one- 
particle states of mass m) we can write the final expression 

/>oo 

-i(0\[(/){x),<j)(x')\0) = ZA(x-x';m 2 ) + / d[i 2 o(p?) A(x - x'; ^ 2 ) (155) 

J m\ 

where the first term represents the one-particle states and the integral represents 
the continuum of multi-particle states with a threshold at mi. In other words, 
if there is a stable particle with mass to, the spectral function must have a 
J-function at to' 2 = to 2 with strength Z, the spectral weight of the one-particle 
state. 

Since the field (f>(x) obeys equal-time canonical commutation relations with 
the momentum IT(a;) = do(f>(x), we get 

-i(0\[n(x,x o ),(j)(x' ,x Q )]\0) = Z lim d A(x - x'; to 2 ) 

r-OO 

+ / dn 2 a(fi 2 ) Jim d A(x - x'\ /j?) 

(156) 

On the other hand, the free field commutator A(x — x';m 2 ) obeys the initial 
condition 



lim d A(x-x';m 2 ) = lim [II (x), 4>{x')] = -iS 3 (x - x") (157) 
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Hence, we find that the spectral function obeys the following identity (known 
as the spectral sum rule) 



/•OO 

l = Z+ dfi 2 a(fj, 2 ) 

J m? 



(158) 



From the positivity of a(m' 2 ) we get, < Z < 1. The lower end of the integra- 
tion range, the threshold for multi-particle production m 2 is equal to 4m 2 since, 
at least, we must create two elementary excitations. A similar analysis can be 




single-particle spectrum 



Figure 15: Spectrum of the propagator. 

done for the Feynman (time-ordered) propagator 

G F {x - x'; m) = -i(0\T(P(x)(p(x')\0) (159) 
which has the spectral representation 

P OO 

G F (x-x';m) = ZG Q (x-x';m)+ dfi 2 <r(n 2 )Go(x - x'; (160) 

This decomposition is known as the Lehmann representation. Here G^(x — 
x';m 2 ) is the Feynman propagator for a free field 

G {x -x') = - / tt-Vj^ 7, (161) 

uv ' J 2tt) 4 p 2 - m 2 +ie K ' 
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In the limit e — > + , the poles of the integrand can be manipulated to give 



lim 



1 



V- 



1 



— in5(p 2 — to 2 ) 



(162) 



e^o+ p 2 — to 2 + ie p 2 — m 2 
Using this identity we get that, in momentum space, the propagator is given by 



G F (p; m) 



p 2 — m 2 + ie 
Its imaginary part is given by 



f-OO 

J mi 



z r . ^(m 2 ) 

p 2 — n 2 + ie 



d^ 2 



Hence 



-Im G F (p; m) = Z 5{p 2 - m 2 ) + a{p 2 )Q{p 2 - to 2 ) 



(163) 



Im G F (p;m) = ttZ S(p 2 - m 2 ) + n a(^ 2 )S(p 2 - ^i 2 ) (164) 



(165) 



pole 




two-particle continuum 



threshold 



cut 



Q 2 



Figure 16: The analytic structure of the propagator is encoded in the spectral 
density a(Q 2 ). 



Once the imaginary part is known, the real part is found through the Kramers 
Kronig or dispersion relation. 



Re Gf(p, m' 



i r c 

k Jo 



dfi 



2 Im G F (p, n 2 ) 
H 2 — p 2 — ie 



(166) 



We see that, in general, there are two contributions to Im G F (p;m). The 
first term is the contribution from the single particle states. In addition, then 
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is a smooth contribution (the second term) which results form multi-particle 
production. While the single-particle states contribute with an isolated pole (a 
5-function in the imaginary part), the multiple particle states (or continuum) 
are represented by a branch cut. 

There is a simple and natural physical interpretation of these results. If the 
incoming state has Q 2 < m 2 it cannot propagate since the allowed value is at 
least, m 2 (the physical mass). If Q 2 > m 2 , and if there are no bound states, the 
incoming state can decay into at least two single particle states. Hence m 2 = 
4m 2 . These states should form a continuum since given the initial momentum 
Pi there are many multi-particle states with the same total momentum. Thus, 
those processes are incoherent. Notice that without interactions, the incoming 
state would not have been able to decay into several single particle states. 

We should stress that the propagators of all the theories that we have dis- 
cussed have the same type of analytic structure that we have discussed here. 



10.5 The S-matrix and the Vacuum Expectation Value of Time 
Ordered Products of Fields 

We are now in position to find the connection between 5-matrix elements and v. 
e. v. of time-ordered fields. For simplicity we will keep in mind the case of scalar 
fields but the results are easily generalizable. The actual derivation is rather 
lengthy and unilluminating. We will discuss its meaning and refer to standard 
texts for details. 

Let's assume that we want to evaluate the matrix element 



<Pi 



,f>„;ou%i, . 



;in) = (pi, . . . ,p n \S\qi, ... ,q m ) 



(167) 



We will assume that all incoming and outgoing momenta are different. This 
matrix element is given by the reduction formula 



(pi,...,p n \S\qi,...,q m ) = 



Z (n+m)/2 



d 4 y x . . . d 4 y n d 4 xx . . . d 4 x n 



x exp 



/ 1 > Pi - ye 

=1 



y^gfc • x k ) 



fc=l 



|2 



k=l 



x(O|T( ( / ) ( 2/1 )... ( / ) (2 / „)0(x 1 )..>(x m ))|O) 

(168) 

where m 2 is the physical mass and the external momenta p and q are on the 
mass shell, p 2 = q 2 = m 2 . 

Let us consider for example the 2 — > 2 process 

(p 1 ,p 2 ;o\it\qi,q 2 ;m) = (p 1 p 2 ,m\S\q 1 ,q 2 ;m) 

= (P1P2, out | <4 n (gi) | <7 2 , in) 

(169) 
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where a\ n (qi) is 



4i(9i) 



-if d 3 x d (j) in {x) e - ^ 1 • x = 

J fixed t 



J d 3 x (uj(q)(j) in (x) + illin(x)) 



e -m ■ x 



(170) 



Hence, the matrix element is 

(pip 2 out|<7i(22hi) = -% lim / d 3 x do{piP20ut\^ in (x)\q 2 m) e~ lqi ' x 

= -i lim — [d 3 x d (p lP 2Out\(t){x)\q 2 m) e~ 1 ^ ' x 



(171) 



where we have made the replacement of by -7^72 4> inside the matrix element 
it -> -00). 
But 

(0|^(x)|l) - ^ 1/2 <O|0 in (^)|l) - ZV^OI^^II) (172) 

and 

( lim - lim ^ f d 3 x F(x, t) = lim lim f ' dt ^- f d 3 x F(x, t) 

\t f ^ + oo U->-ooJ J tf^+OOti^-OO J t , dt J 

(173) 

These formulas allow us to write 

lim / d 3 x (pip 2 ,out|aJ n (qi)|g 2 ,in) = 

tf -^+oc J 

= ^\m\^jd 3 x (pip 2 ,out|a[ n ((7i)|(72,in) + J d A x d (piP2, out 1 0^(91)^2, in) 
= (PiP2,out\al ut (qi)\q 2 ,'m) 

(174) 

Thus, the matrix element is 

(piP2,out|<7i<72,in) = (piP2,out|4 ut (gi)|q 2 in) 



+ 



Z 1 / 2 



J d 4 x do 



e lqi ' x d {pip 2 out\<j>{x)\q2m) 



(175) 



The first contribution is a disconnected term and it is given by 
(piP20ut|a£ ut (gi)|g 2 in) 

= {2n) 3 2p° 1 5 3 (p 1 - gi)<P20ut|giin) + {2tt) 3 2 P 2 5 3 (p 2 - gi)<Piout|giin) 

(176) 



34 



Notice that, if qi is on the mass shell, q\ — m 2 , and e l 1 lX is a solution of the 
Klein- Gordon equation 



(d 2 + mV?i ' x = (qf^m 2 ) 



(177) 



The second contribution to the matrix element can also be evaluated (for arbi- 
trary states a and (5 ) 



d 4 xd Q e- lqi - x do{p,ont\4>{x)\a,in) 
= J d 4 x [e - ^ 1 ,a; ^(/3,out|0(a;)|a,in) + (-^e - ^ 1 ' x )(f3 1 out|0(x)|a, in) 
= y d 4 x e- i<zl - x ^(/3,out|^(a;)|a,in) + (- V 2 +m 2 ) e - ^ 1 • » ut|^(x)|a, in) 
= J d 4 xe- 1 ^ ■ x (d 2 + m 2 )(/?,out|<£(a;)|a,m) 



(178) 



where we have integrated by parts. Hence, the matrix element is 



(f>ip 2 ,out|gig 2 ,ii 



Z 1 / 2 



d x\e 



(d +m )(pip2,out\<j}(xi)\q 2 ,m) 



+ {2ir) 6 2p\d i {p 1 -q 1 ) (p2,out|ga,in) + {2Trf2p 2 > 8'\p 2 - «i)<Pi, ou% 2 , in) 

The matrix element inside the integrand of the first term is equal to 

(pip 2 ,out|0(xi)|<? 2 ,in) = lim — ^ / d 3 yie tpi ' Vl d y o(p 2 ont\(f>(y 1 )\q2m} 

j/J->+oo Z ' J 

(180) 

where (by definition) yf > x\. This expression is also equal to 
(PiP20ut|</>(xi)|g 2 in) = (p 2 ,out|0(x 1 )a in (pi)|q 2 ,in) 



(179) 



+ 



ZV2 



d%e^ ■ y^(d 2 yi +m 2 )(p 2 ovit\Tcb(y 1 )cf,(x 1 )\q 2 m) 



(181) 



By substituting back into the expression from the matrix element we find that 
the latter is equal to 

(pip 2 ,out|<7i<7 2 ,in) = 

(2tt) 3 2p°5 3 ( Pi - qi )(p 2 out\q 2 m) + (27r) 3 2p^ 3 (p 2 - gi)(piOut|g 2 in) 



+ 



Z 1 / 2 



d A Xl e~ 1 ^ ■ x i (d 2 Xl + m 2 )(p 2 out|0(a ;i |0, in) {2^f2q° 2 5 3 {q 2 - Pl ) 



+ (^2) 2 j d 4 xid* yi e^ ■ Vi ~ 91 • *i) (d 2 Xi + m 2 ){d 2 yi + m 2 ) 
x (p2out|r0(2/i)0(xi)|g 2 in) 



(182) 



35 



By iterating this process once more we obtain the reduction formula plus dis- 
connected terms. The reduction formula provides the connection between the 
on-shell S'-matrix elements and v. e. v. of time ordered products. Notice that 
the reduction formula implies that the Green functions must have poles in the 
variables p\ (where pi is conjugate to X{) and that the matrix element is the 
residue of this pole. We will see later on that this residue is the on-shell one 
particle irreducible vertex function. 

The reduction formula shows that all scattering data can be understood in 
terms of an appropriate v. e. v. of a time ordered product of field operators. 
The problem that we are left to solve is the computation of these v. e. v. 's. We 
will use perturbation theory to find them. 

10.6 Linear Response Theory 

In addition to the problem of evaluating S'-matrix elements, it is of interest 
to consider the response of a system to weak localized external perturbations. 
These responses will tell us much about the nature of both the ground state 
and of the low-lying states of the system. This method is of great importance 
for the study of systems in condensed matter physics. 

Let H be the full Hamiltonian of a system. We will consider the coupling 
of the system to weak external sources. Let A(x,t) be a local observable such 
as the charge density, the current density or the local magnetic moment. Let 
us represent the coupling to the external source by an extra term H ext (t) in the 
Hamiltonian. The total Hamiltonian now is 

H T = H + H cxt (183) 

If the source is adiabatically switched on and off, then the Heisenberg represen- 
tation for the isolated system becomes the interaction representation for the full 
system. Hence, exactly as in the interaction representation, all the observables 
obey the Heisenberg equations of motion of the system in the absence of the 
external source while the states will follow the external source in their evolution. 
Namely, let |gnd) be the exact ground state (or vacuum) of the system in the 
absence of any external sources. The external sources perturb this ground state 
and cause the v.e.v. of the local observable A(x,t) to change: 

(gnd|i(aU)|gnd) -> (gndlC/" 1 ^) A(x,t) U(t)\gnd) (184) 

where the time evolution operator U(t) is now given by 

i /"* 

U{t)=Texp{--J dt'H ext {t')} (185) 

Linear Response Theory consists in evaluating the changes in the expectation 
values of the observables to leading order in the external perturbation. Thus, to 
leading order in the external sources, the change is 

i r* 

6(gnd\A(x,t)\gnd) = - dt' (gnd\[H ext (t'),A(x,t)]\gnd) + ... (186) 
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Quite generally, we will be interested in the case in which H oxt represents the 
local coupling of the system to an external source f(x, t) through the observable 
A(x,t) . Thus, we will choose the perturbation H ext (t) to have the form 

tfext(i) = J d 3 x f(x,t) A(x,t) (187) 

The function f(x,t) is usually called the force. If the observable is normal or- 
dered relative to the ground state of the isolated system, i.e., (gnd\A(x, t)|gnd) = 
0, the change of its expectation value will be equal to the final value and it is 
given by 

6(gnd\A(x,t)\gnd) f = l -j dt' J d 3 x' (gnd\A(x,t'),A(x,t)]\^d) f(x',t') + --- 

(188) 

The main assumption of Linear Response theory is that the response is propor- 
tional to the force. The proportionality constant is interpreted as a generalized 
susceptibility \- 

<5<gnd|i(f,i)|gnd> = X -f= f d 3 x' f dt 1 X (x,x') f(x') + ... (189) 

where x( x i x ') 1& the susceptibility. 

Let D R (x, x') represent the retarded Green function of the observable A(x, t), 

D R (x,x') = -ie(ar -a:J,)<gnd|[i(a:),i(a;')]|gnd) (190) 

We see that (A(x)) is determined by D R (x,x') since 

(5(gnd|^(a;)|gnd) = i J d 4 x'D R (x, x')f(x') + . . . (191) 

Therefore, the responses and the susceptibilities are given by retarded Green 
functions, not by time-ordered ones. However, since the retarded and time- 
ordered Green functions are related by an analytic continuation, the knowledge 
of the latter gives the information about the former. 

Let us Fourier transform the time-dependence of (A(x,t)). The Fourier 
transform, (A(x, u>)}, is given by the expression 

5(gnd\A(x,w)\gnd) = 

J d*x'{- % - £ dr(gnd\[A(x,t),A(x' ,t + r)]|gnd) e wr |/(x',c) 

(192) 

where f(x,oj) is the Fourier transform of f(x,t). Thus, the Fourier transform 
of the generalized susceptibility, x( x ,x';uj) is given by 

X (x,x" ;uj) = ~ [ dre lLUT (gnd\[A(x,0)A(x" ,r)]|gnd) (193) 
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which is known as the Kubo Formula. Hence 



I 



— OO 



dre llVT D R (x,x') 



(194) 



If we also Fourier transform the space dependence, we get 



(A(p,uj)) = \D R {p,uj)f{p,uj) 



(195) 



and the generalized susceptibility x(p, uj) now becomes 




(196) 



In practice we will compute the time- ordered Green Function D(x,x'). If we 
recall our discussion about the propagator, we expect D(p, ui) to have poles on 
the real frequency axis. For D(x,x') to be time ordered, all the poles with 
uj < should be moved (infinitesimally) upwards into the complex frequency 
upper half-plane while all poles with uj > should be moved downwards into 
lower half-plane. Thus D(p, u>) is not analytic on either half-plane. The retarded 
Green function D R (p, uj) is (with our conventions for Fourier transform) analytic 
in the lower half-plane. Thus we can relate D(p,ui) to D R (p,uj) by 



The time-ordered Green function D{p,ui) (the propagator for the observable 
A(x,t)) admits a spectral (or Lehmann) representation similar to that of the 
propagator for the relativistic scalar field, Eq.(160). Similarly, we can define the 
spectral function A(p, uj) of the observable to be A{p,ui) = ImD rot (p,w). 

These relations imply that the susceptibility x(p> w ) obeys the Kramers- 
Kronig (or dispersion) relation 



Finally, let us recast the formulas for a general change of an arbitrary operator 
into a more compact form. We can apply the formulas that we derived for the 
interaction representation just to the part of the Hamiltonian which involves 
the coupling to the external sources ff ext (i). The interaction representation 
S'-matrix is 



ReD R (p, uj) = ReD(p, uj) = hRex(p, u>) 
lmD R (p,ui) = lmD(p, uj) signuj = hlmx(p,u)) 



(197) 




(198) 




(199) 




gnd) (200) 
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which we denote by Z[f] 



Z[f] = (gnd|T e 




(201) 



By expanding in powers of f(x) we get 

Z[f] = 1 - l - J d 4 xf(x){gnd\A(x)\gnd) 

+ h { ~Ji )2 J d " X J dix 'f( x )f( x ')^d\TA(x)A(x')\gnd) + ... 

(202) 

The second term vanishes if A(x, t) is normal ordered. Within the same degree 
of precision, we can re-exponentiate this expression to give 

^zfd'xf d 4 x'f(x)K(x, x')f(x') + 0(f) 
Z[f] = e 2h J J (203) 

where the kernel K (x, x') is the time-ordered Green function 

K(x,x') = Ugnd\TA(x)A(x')\gnd) (204) 
n 

In general the observables of physical interest are, at least, bilinears in the fields. 
Thus, the kernels K(x,x') represent not one-particle propagators but generally 
propagators for two or more excitations. We can learn a lot from a physical 
system if the spectral functions of the kernels K(x,x r ) are known. In general 
we expect that the spectral function will have a structure similar to that of 
the propagator: one (or more) <5-function contributions and a branch cut. The 
J-functions are two (or more) particle bound states known as collective modes. 
The branch cuts originate from the two or multi-particle continuum. Examples 
of collective modes are plasmons (sound waves) in electron liquids and spin 
waves in magnets. 



10.6.1 Application of the Kubo Formula: Electrical Conductivity of 
a Metal 

We will only consider the response of an electron gas to weak external electro- 
magnetic fields A^(x). The formalism can be generalized easily to other systems 
and responses. In particular, we will consider the electrical conductivity of a 
metal. 

There are three effects (and couplings) that we need to take into consider- 
ation: a) electrostatic, b) diamagnetic (or orbital) and c) paramagnetic. The 
electrostatic coupling is simply the coupling to an external potential with H eyA 
given by 

#ext= Yl [ d 3 x e<f>(x,t) i>l(x)i>„{x) (205) 

c=t>i 
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where <f> = Aq is the scalar potential ( or time component of the vector potential 
A^). The diamagnetic coupling (or orbital) follows from the minimal coupling 
to the external vector potential A. The kinetic energy term Hk is modified 
following the minimal coupling prescription to become 



which can be written as a sum of two terms 



x))ipUxy( V --A(x))Mx) (206) 



H K (A) = H K (0) + H eKt {A) 



(207) 



where Hk(0) is the Hamiltonian in the absence of the field and H ext (A) is the 
perturbation, i.e., 



H cxt (A) = J d 3 x J{x) ■ A(x) - ^A 2 (x) £ V£(*)Vv(z) 



(208) 



Here J(x) is the gauge-invariant charge current 

Zmc z — ' L J mc z * — ' 



ieh 
2mc 



ipt(x)DMx)- (D^{x)) ^{x) 



(209) 



Clearly J(x) is the sum of the two terms, one which represents the mass current 
and the diamagnetic term, A 2 J2 a V'JVv- We can write the total perturba- 
tion, including the scalar potential A , if we write 



d 3 G 



Mx)A»(x) 



2mc 2 



(210) 



Finally, we can also consider a paramagnetic coupling to the spin degrees of 
freedom which has the Zeeman form 



ti, 



Zccman 



f d 3 xgB(x)-J2^)Sa^a'(x) (211) 



where g is typically of the order of the Bohr magneton [Lb and S = |<r for spin 
one-half systems. 

A straightforward application of the Linear Response formulas derived above 
yields an expression for the current (J^)' in the presence of the perturbation. 



dt' / d 3 x'(gnd\[J„(x'),J^x)}\gnd)A,(x') (212) 
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This formula suggests that we should define the retarded current correlation 
function T>j} l/ (x, x') 

V^ t (x,x , ) = -iQ(x -x , )(gnd\[J l ,(x),J l/ (x , )]\gnd) (213) 
The induced current (J M )i nc j 

Wind = W - (j^gnd (214) 

(where j M is the mass current) has a very simple form in terms of T>j} u (x,x'), 
namely 

(Mx)) ind = l f d 4 x'V*?(x,x')A»(x')- -ilA fc (aO<n(aO)<W +0(A 2 ) (215) 

Below we will show that (J M (x))i n d is (a) conserved (i.e., d*(J^(x))i n< i = 0) and 
(b) gauge-invariant. 

Since (J Ai (x)); n d is gauge invariant, we can compute its form in any gauge. 
In the gauge Ao = the spatial components of (J M (x))i n d are 

(J fc (x)) ind = -^A k {x) + / d 4 x'Vlf(x - x')A e (x') + 0(A 2 ) (216) 



In this gauge, the external electric field E ext and magnetic field H are 

E ext = -8 A H = \}xA (217) 
Now, in Fourier space, we can write 

2 

(J fe (p») ind = -^A k (p,u) +V I k f(p,Lu)A e (p,Lu) 



2 \ T^cxt 

7>$tf,u)--%r5 kt )-?-(p,u) 



iuj 



(218) 



This expression is almost the conductivity. It is not quite that since the con- 
ductivity is a relation between the total current J = J; nd + J C xt and the total 
electric field E. In order to take these electromagnetic effects into account, we 
must use Maxwell's equations in a medium, which involve E, D, B and H 



S7-D = p vx£ = -f 



dE 

V • J-> — u V^-"— ~~ 

where 



(219) 



B = H + M E = £ cxt + E ind (220) 

Here M and E lnd are the magnetic and electric polarization vectors. In partic- 
ular 

J md = d t E iad (221) 
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and 

d t D = 8 t E + J md (222) 
Linear Response theory is the statement that D must be proportional to E, 

Dj = EjkEk (223) 

where Sjk is the dielectric tensor. Since E and E cxt satisfy similar equations 

-V x V x E = d 2 E + d t J 
-V x V x # cxt = d 2 E cxt + d t J cxt 

(224) 

Since y x y x E = y(y ' E) — X/ 2 E, we can write, for the Fourier transforms, 
the equations 

PiPjEj(p, uj) - p 2 E l {p, w) = -Lu 2 Ej(p, w) - iLuJi(p, w) 
PlPj E] xt (p, w) - fEl xt (p, u) = -uj 2 E^ (p, w) - iu Jf xt (p, uj) 

(225) 

Thus, we get 

PiPjEj(p, w) - p 2 E l (p, uj) + uj 2 Ei{p, uj) = -iujJl nd (p, w) 

+ p iPj Ef (p, uj) - p 2 E?* (p, uj) + uj 2 Er l (P, lo) 

(226) 

and 

-iojJ? d (p, uj) = [Sij^ - V«(p, u)j E^(p, uj) (227) 
from where we conclude that 

(PiPj - P 2 5ij + ) Ej (p, uj) = 

2 

{Si^ - V«j(p, uj) + PiPj - p 2 5 tJ + uj 2 8 10 ) Ef(p, uj) 

(228) 

In matrix form, these equations have the simpler form 

(p ® p - fl + uj 2 I)E = ( -^4 J - 2? fl + p ® p - p 2 7 + c^ 2 /) £ cx t (229) 

This equation allows us to write -E ext in terms of E. We find that the induced 
current is 

iujJind = 

{p®p-p 2 I + uj 2 I) E 
(230) 



V 1 



e p 
mc 2 



rac^ 



■p (X p — p ' 



■to 2 1 



42 



and we find that the conductivity tensor a is 



e 2 p 
mc 2 



I + 



9 



mc' 



2 



— %J -V H (p,u)+p(g>p-p 2 I + lo 2 I 



mc' 



Also, since £> = eE, the dielectric tensor e is 



(231) 



e = 7+ - <t 



(232) 



10.6.2 Correlation Functions and Conservation Laws 

In the problem discussed in the previous section, we saw that we had to consider 
a correlation function of currents. Since the currents are conserved, i.e., 



(233) 



we expect that the correlation function T>^ v (x,x') should obey a similar equa- 
tion. Let us compute the divergence of the retarded correlation function, 



d;vZ(x,x') = dl (-ie(x - a/ ){gpd\[J»{x), JV)]|gnd» 



(234) 



Except for the contribution coming from the step function, we see that we can 
operate with the derivative inside the expectation value to get 



&P%{x,a>) = -i(d^e(x Q -x' )) (gnd\[J»(x),J»(x')}\gnd) 
-ie(x -x' ){gnd\[d^(x),r(x')}\gnd) 



(235) 



The second term vanishes since J M (x) is a conserved current and the first term 
is non zero only if (i = 0. Hence we find 

&PZ{x,x') = -i6(x - 4)(gnd| [J°(x), r(x')] |gnd) (236) 

which is the v.e.v. of an equal-time commutator. These commutators are given 
by 

(gnd| [J°(x,x ),J°(x,Xo)] |gnd) = 



le 



(gnd| [J°(x,x ),J i (x',xo)] |gnd) = ^d%[5{x - f)<n(f))] 



Hence, the divergence of V 1 ^ is 



(237) 



d»V$(x,x') 



mc- 



;d* k [S\x-x')(n(x))] ■ d^(x,x') = (238) 
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and 

d v x ,V%{x,x') = -^d%[5 i {x-x'){n{x'))] ■ d x i ,V$(x,x')=0 (239) 

Notice that the time-ordered functions also satisfy these identities. These iden- 
tities can be used to prove that ( J md ) is indeed gauge-invariant and conserved. 
Furthermore, in momentum and frequency space, the identities become 

-iwX>oo (p, w) - ipkDka (P) w ) = 
-iwV§ (p, u) - ip k T>o k (p, u) = 



9 



-iuT>k (p, w) - ipiD^ (p» = -ip k 



(240) 



We can combine these identities to get 



uj 2 vH(p, u) - pip k Vf k (p, u) = -^p 2 (241) 

Hence, the density-density and the current-current correlation functions are not 
independent. A number of interesting identities follow from this equation. In 
particular if we take the static limit u) — > at fixed momentum p, we get 

2 - 

lim impkVlf (p, u) = e -^p 2 (242) 

provided that lim w ^ 2? oo*® w ) i s not singular for p ^ 0. Also from the equal- 
time commutator 

(gndp^xoMo^a^lgnd) = ^8% (S(x - x' )(n(x))) (243) 



we get 



km d%Vlt{x,x') = ^ V l{5{x-x'){n{x))) (244) 



If the system is uniform, (n(x)) = n = p, we can Fourier transform this identity 
to get 

r+oo , 

ip k V k t(p,u>) = -^p 2 (245) 



The conservation laws yield the alternative expression 



/ 



L 



2ir mc 
which is known as the f '-sum rule. 



<! "iu;V$(p,u) = ^p 2 (246) 
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If the system is isotropic, these relations can be used to yield a simpler form 
for the conductivity tensor. Indeed, for an isotropic system T>j$(p,uj) can only 
have the form of a sum of a longitudinal part Vl ct and a transverse part 2?™* 



Thus, we get a relation between and 2?jj ct 



V P 2 



(247) 



e n ^ 9 
2P 



Hence 



and 



2?f(p, W ) 



2 - 
en 



2 - 
en 



(248) 
(249) 
(250) 



liml?r(p» = 

w^o 11 v mc 2 

for all p. 

The conductivity tensor can also be separated into longitudinal o\\ and trans- 
verse cr_L pieces 



u - °ll -~T + IFF ~ dl 



P 



P1P3 



we find 



and 



->rct 



e 2 n 



q~yrct e n 

II mc 2 



_£>rct 



1 



mc 2 



1~yret en 

-L mc 2 " 



mc 2 



■ V r f +LU 2 -p 2 



(251) 
(252) 

(253) 



These relations tell us that the real part of <tii is determined by the imaginary 
part of ^jj *- Thus, the resistive part of <j\\ (which is responsible for dissipation 
in the system) is determined by the imaginary part of a response function. This 
is generally the case. 



A The Dirac Propagator in a Background Elec- 
tromagnetic Field and Coulomb Scattering 

Let us consider briefly the Dirac propagator in a background electromagnetic 
field and use it to compute the S'-matrix for Coulomb scattering. By a back- 
ground field we mean a classical (fixed but possibly time-dependent) electro- 
magnetic field A^{x). We will denote by Sf{x,x'\A) the Dirac propagator in a 
background field A^. Sf{x,x'\A) obeys the Green function equation 

{i$ -eA-m) S F (x, x'\A) = S 4 (x - x') (254) 
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whereas the Dirac propagator Sf(x,x') in the absence of a background field 
obeys instead 

m) S F (x,x') = S 4 (x- x') (255) 
Thus, we can also write Eq.(254) as 

(Sp 1 -e£)S F (A) = l (256) 

Hence 

S F (x, x'\A) =S F (x-x') + eJ d 4 y S F (x - y) A{y) S F {y, x'\A) (257) 
or, in components, 

Sf(x,x'\A) = Sf(x-x') + e J d 4 yS a F \x-y) [A^y)^ Sf{y,x'\A) 

(258) 

As an explicit application we will consider the case of Coulomb scattering of 
Dirac electrons from a fixed nucleus with positive electric charge Ze. We will 
now compute the S'-matrix for this problem in the Born approximation. As in 
non-relativistic Quantum Mechanics, in this approximation we replace the prop- 
agator in the integrand of Eqs. (257) and (258) by the free Dirac propagator, 
S F (x-x'). 

Consider now an incoming state, a spinor the we will denote by ^i(x), with 
a particle with positive energy (an electron) and spin up (say in the z direction), 
and momentum pi. This incoming (initial) state is (for x — > — oo) 

*i(x) = u< a >(pi) ^e-'Pi-x (259) 



V 

The outgoing (final) state / (x) is a spinor representing also a particle with pos- 
itive energy (an electron) with spin up (also in the z direction) and momentum 
Pf, and it is given by 

* J (y) = l u W( pf )[™l e -iPf-V (260) 



w 

The S'-matrix is 

S fi =i lim lim d 3 x d 3 y^ f {y,y ) S F (y,x\A) Vi(x,x ) (261) 

xo— > — oo j/o— >+oo J J 

At the level of the Born approximation we can write 

Sn = i lim lim d 3 x d 3 y^ f (y,y ) S F (y,x) ^i(x,x ) + 

lim lim d 3 x d 3 y d 4 z^ f {y 7 y )S F {y 7 z)Mz)S F (z 7 x)^ l {x,x ) + ... 

o^-ooyo^+oc J J J 



I 

xo- 



(262) 
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We now recall the expression for the free Dirac propagator 

S F (x -x' = -i(0\Ti> a (x)i> a ,(x')\0) 
d 3 p 



(2 



( e ( x, o- x o)e- i P^ x '-^A + (p) + 9(x -x' )e- i P-( x - x 'A.(p)) 

(263) 



where A±(p) are projection operators onto positive (particle) and negative (anti- 
particle) energy states: 

A ± (p) = -J-(±jrf + m) (264) 
Ira 

Alternatively, we can express the propagator in terms of the basis spinors u a (p) 
(which span the positive energy states), and v a (p) (which span the negative 
energy states), as 

S F (x'-x)= -iO(x' -x ) I d 3 p ^\x')u^\x) 

J (7=1,2 

+i6(xo-x' ) f d 3 p £ v<f\x')v^(x) 

J (7=1,2 

(265) 

where we have used the notation 

u^\x) = u^\p)e- i P- x , v'f\x)^v^{p)e- % P- x (266) 

Let us begin by computing first the top line in Eq.(262), the projection of the 
free propagator onto the initial and final states. By expanding the propagator 
we find, 

j d 3 xd s yi> f {y)S F {y-x)i) i {x) = 

J (7=1,2 

+i6(x -y ) J d 3 p J d 3 xd 3 yj f (y)v^(y)v^(x)Mx) 

(j— 1,2 

(267) 

We now use the orthogonality relations of the Dirac basis spinors to find 
J d 3 y^f\y)u^(y)=S^6 3 (p-p f ) 
J d 3 x u^\x)4 a \x) = 8° a 8 3 {p-p l ) 

(268) 
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Hence, to leading order the matrix element S/j of the S'-matrix is 

S fi = 5 3 {p f - pi)5 af) + Born term (269) 

Let us now compute the Born term (the first Born approximation). We will 
need to compute first an expression for 

J d 3 yj f (y)S F (y,z) (270) 

and for 



d 3 x S F (z,x)il>i{x) (271) 
Using once again the expansion of the propagator we find that Eq.(270) is 

d 3 y iP f {y)S F {y,z) = 



/' 

J d 3 y^f\y)(-i)9(y -z ) £ J d 3 p «W (y)u^(z) 

(7— 1,2 
(7 = 1,2 

ie( yo -z Q )j2 I d3 p ( I d3 y $f (»)<(»)) 



/ 



+i0(z„ - 2/ ) ^ fd 3 p f f d 3 y $f(y)v?(v)) t\ z ) 

{—i9(yo — zo) rpf(z), if the final state is a particle 
+i9(z a — y ) ipf(z), if the final state is an antiparticle 



(272) 



The other expression, Eq.(271), can be computed similarly. Putting it all to- 
gether we find that the Born term is 

Born term = -ie J d A z i>f\z) A^\ a \z) %o - zo) d{z - x ) (273) 

corresponding to an electron propagating forward in time. 

Let us evaluate this expression for the case of a Coulomb potential, 

A» = (A ,0), A = ^ (274) 

with r = \z\. The Born term now becomes 

Born term = -ie J dz Q J d 3 z $p 1 (z^otp^ ' (z) 9(y - z ) {^^j 

= ^ r dz [ d 3 z e <(P/ "ft) • * 1 uW( Pf ) l0 u^( Pl ) 

V y/EiEf 47T J_ OQ J r 

(275) 
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where V is the volume. Using now that 

r dz Q e^ E f ~ E ^ z ° = 2ir5(E f - E,) (276) 

J — oo 

we can write the Born term as 

Born term = — ^_ — _ - J d»r ^>( P/ )W^) 

(277) 

2 

where a = |^ is the fine structure constant, q = pf — pi is the momentum 
transfer, and 

, e- i $- r 4tt 
r \q\ 

The matrix element of the S'-matrix, in the Born approximation, is then 
equal to 

S fi = 6^6 3 (p f -p l ) + l ^ -JL= 2*5(E i -E f ) % «(«(p / )7o« (a) (p i ) (279) 
Since 

# states withp) within d 3 p f = V^-^r (280) 

(27T) J 

we can write the transition probability per particle into these final states as 

lSfil v w = w { f l) 

(281) 

where T is the time of measurement (this is Fermi's Golden rule). 
Thus, the number of transitions per particle and unit time is 

- / if im, w -"w-w - < 282 ) 

Dividing out this expression by the incoming flux, y ^ , we obtain an expression 
for the differential cross section 

da fi = (| dp f p) ^^\u^( Pf h u^( Pl )\H(E f - E^dQf (283) 

For elastic scattering, \pi\ — \pj\ — pj and EdE = Pfdpf, we obtain that the 
differential cross section is 

4Z 2 a 2 m 2 
For an unpolarized beam we get 



*>ti = ^ |^ W (p/)7o^ (a) (p»)l 2 ^/ (284) 



where j3 — v/c. 



Z 2 a 2 



unpolarized 4\p\ 2 ft 2 Sm 4 (0/2) \ 2 



1 - /3 2 sin^ - (285) 
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